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ABSTRACT 

We report FUSE 1 and new deep optical spectroscopic observations of PNG 135.9+55.9, the 
most oxygen-poor planetary nebula located in the Galactic halo. These observations allow us 
to estimate the gravity of the central star by fitting the profile of the observed H I absorption 
lines with NLTE model atmospheres. Our best fit implies that the central star is still in a 
pre-white dwarf stage. We also find large variability of the radial velocities of the absorption 
component of the Balmer lines on a timescale of hours. This is direct evidence that the nucleus 
of PNG 135.9+55.9 is a close binary. The large semi-amplitude of the radial velocity variations 
and the probably short period suggest a massive compact companion, likely a white dwarf. 
Although our orbital solutions are very preliminary, they indicate that the total mass of the 
system probably exceeds the Chandrasckhar limit. If confirmed, this would make this binary the 
potential progenitor of a type la supernova. 

Snhipst h.p.n.dinns: TSM: nla.neta.rv nfihulae: individual ofmt.ra.l star: Ultraviolet: TSM: snfictrnsmnv: 



1. Introduction 

PNG 135.9+55.9 (also called SBS 1150+599 A) 
is a recently discovered planetary nebula in the 
Galactic halo (Tovmassian et al. 2001) whose oxy- 
gen abundance is reported to be extremely low, 
of order 1/100 of the solar value or less (Tov- 
massian et al. 2001; Richer et al. 2002; Jacoby et 
al. 2002). This makes it, by one order of mag- 
nitude, the most oxygen-poor planetary nebula 
known so far. Such a low oxygen abundance may 
be either genuine and represent the oxygen abun- 
dance in the material from which the progenitor 
star formed, or it may be the result of mixing pro- 
cesses while the star was on the AGB. Although 
the presently available determinations of the Ne/O 
ratio are discrepant by a factor of about 10 (Richer 
et al. 2002 and Jacoby et al. 2002), they indicate 
that, even in the case of conversion of O into Ne, 
the original oxygen abundance of the progenitor 
of PNG 135.9+55.9 is likely extremely low. 

Clearly, this unique object deserves a more de- 
tailed study to better define its nature. Tovmas- 
sian et al. (2001) speculated upon various scenar- 
ios for the formation of this planetary nebula. One 
of them assumed a progenitor star having formed 
out of infalling material. Another invoked the ex- 
istence of a close binary core. In this paper we 
focus on the information provided by new obser- 
vations concerning the nature of the central star of 
PN G 135.9+55.9. A companion paper (Stasihska 
et al. in preparation) will discuss new observa- 
tional constraints on the chemical composition of 
this object. 

First, we acquired far- ultraviolet spectra of 
PNG 135.9+55.9 using the Far Ultraviolet Spec- 
troscopic Explorer (FUSE). We also obtained a 
series of high signal-to-noise optical spectra of the 
nebula and its central star. From these obser- 
vations, one can constrain the fundamental pa- 
rameters of the central star. It is also possible 
to investigate whether there might be variability 
on short time scales. The most surprising re- 
sult of our study is the discovery that the core of 
this planetary nebula is a short period, close bi- 
nary with an unseen companion, probably a white 
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dwarf. 

The paper is organized as follows: In Section 
2, we describe the observations and the reduc- 
tion procedure, both for the FUSE and ground- 
based observations. In Section 3, we describe 
the main spectral features of the object from the 
FUSE and ground-based observations. In Section 
4, we present a stellar atmosphere analysis that 
allows us to infer the gravity of the ionizing star 
from the observed H I absorption features and es- 
timate parameters of the binary system. In Sec- 
tion 5, we present estimates of the reddening of 
PNG 135.9+55.9 based upon the FUSE data and 
in Section 6 we estimate the distance and age of 
PNG 135.9+55.9. In Section 7, we discuss the 
binary nature of the core and the masses of its 
components. Finally, in Section 8, we summarize 
the main results of this study and present some 
prospects for future study. 

2. Observations 

2.1. FUSE observations 

FUSE is an orbiting observatory with the ca- 
pability to obtain high resolution spectroscopy in 
the wavelength region between 900 and 1200 A. It 
consists of four co-aligned telescopes optimized for 
far-UV wavelengths. The spherical, aberration- 
corrected, holographic diffraction gratings dis- 
perse light from four channels. Two channels with 
SiC coatings cover the 905- 11 00 A region and two 
others with Al+LiF coatings detect the 1000- 
1187A range. An overview of the FUSE mission 
has been given by Moos (2000). 

The observation of PN G 135.9+55.9 (ID# 
C034) was conducted through the LWRS aper- 
ture (30" x 30" square) on 2002 January 30 for 
a total exposure time of 30 660 s. The observa- 
tion was composed of seven individual integrations 
of approximately 4200 s each, simultaneous in all 
channels. These observations were among the first 
executed after the recovery from the failure of two 
reaction wheels in December 2001. Therefore, 
the pointing was not as good as was achieved be- 
fore the accident or after the situation was taken 
under complete control. We have checked pos- 
sible deviations during the observations. There 
was a systematic, quasi-sinusoidal pointing drift 
on an orbital time-scale. However, the maximum 
displacement is only 2" 3, which corresponds to 
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a maximum shift of 0.026 A (i.e., 7.7 km/s at 
1000 A), or about half a spectral resolution ele- 
ment. Thus, it can hardly pose a problem, but the 
pointing itself (centering in the aperture) , which is 
difficult to verify, can introduce zero-point offsets 
of up to 0.15A. The pipeline reduction package 
CalFUSE 2.1.7 was used for the extraction of the 
spectra as well as reduction and calibration pur- 
poses. 

Our examination of the spectra showed that the 
wavelength shifts between individual exposures, 
which are common for FUSE data, never exceeded 
1.5 pixels and were usually in the 0.1-0.5 pixel 
range. Since the shifts were negligible, we used 
a non-standard processing wherein all 7 exposures 
were concatenated into one large photon list before 
pushing it through CalFUSE 2.1.7. This permits 
a more reliable background model to be created, 
which is especially important for faint targets such 
as PNG 135.9+55.9, and also eliminates the need 
to cross-correlate and coadd exposures. The other 
important procedure that we used was to separate 
the night exposure (NE) from the total exposure 
of the object. The data are severely contaminated 
by the airglow lines arising form the illuminated 
atmosphere of the Earth during day time (Fcld- 
man et al. 2001). In order to eliminate the strong 
emission lines from the terrestrial airglow, we ex- 
tracted photons obtained when the exposure was 
taken during orbital night. This decreased our ex- 
posure time by factor of 3. We use both the total 
exposure and NE spectra to analyze the data of 
PN G 135.9+55.9. No event bursts or other unusu- 
ally strong glitches common for FUSE data were 
detected in the raw images or the photon arrival 
rate plots (see The FUSE Instrument and Data 
Handbook edited by Sahnow et al. 2002). 

2.2. CFHT optical observations 

The observations at the Canada-France-Hawaii 
Telescope (CFHT) were obtained on 2003 May 
1 UT with the MOS spectrograph (Le Fevre et 
al. 1994). The U900 grism was used to ob- 
serve the 3400-5300 A spectral interval. A 1" 
slit was used. Coupled with the EEV1 CCD, 
2048 x 3900 x 13.5/im pixels, the dispersion was 
0.78A/pix and the resolution 3.0-3.5 A, based 
upon the widths of arc lamp lines. A B filter 
(CFHT filter #1412) was used to center the cen- 
tral star of PN G 135.9+55.9 in the slit. The spec- 



tra of PN G 135.9+55.9 were obtained in pairs with 
the spectrograph slit set at the parallactic angle 
corresponding to the midpoint of the pair of spec- 
tra. This was done since the spectrograph rotation 
and object re-acquisition took a significant amount 
of time. For the observations of PN G 135.9+55.9 
and the standard star BD+33°2642 (Oke 1990), 
the slit position angle was within 10° of the par- 
allactic angle at all times. Seven 1800 s expo- 
sures of PNG 135.9+55.9 were obtained spanning 
the airmass range 1.30-2.65. The observations of 
BD+33°2642 were obtained using a 5" slit. Spec- 
tra of Hg, Ne, and Ar lamps were obtained for 
wavelength calibration at the end of the night only. 
Images of the flat field lamp were obtained through 
both slits to correct for pixel-to-pixel variations. 
Bias images were also obtained in order to remove 
any two-dimensional structure in the bias level. 

The data were reduced using the Image Reduc- 
tion and Analysis Facility (IRAF) 2 software pack- 
age (specifically the specred package). The mean 
level of the overscan region was subtracted from 
each image. Next, the average of the overscan- 
corrected bias images was subtracted from all ob- 
ject and flat field images to remove any two- 
dimensional bias structure. The flat field images 
were combined and then processed to divide out 
the shape of the flat field lamp. The resulting flat 
field was then divided into all object images. The 
spectra of PNG 135.9+55.9 and the standard star 
were extracted to one-dimensional images, sub- 
tracting the sky contribution by defining sky aper- 
tures on both sides of the object spectra and inter- 
polating between them with a straight line. The 
wavelength axis was calibrated using the spectra 
of the arc lamp. Since the object was observed 
repeatedly over a large airmass interval, an em- 
pirical atmospheric extinction curve was derived 
from the data. This extinction curve turned out 
to be somewhat steeper than that given in the 
CFHT Observatory Manual 3 . The spectra were 
calibrated in flux using the observations of the 
standard star and the previous empirical atmo- 
spheric extinction curve. The individual spectra 
of PN G 135.9+55.9 were completely calibrated be- 

2 IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Associated Uni- 
versities for Research in Astronomy, Inc., under contract 
to the National Science Foundation. 

3 See http:/ /www. cfht.hawaii.edu/Instrumcnts/ObscrvatoryManual/ 
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fore combining or further analyzing them. 

3. The main features of the observed spec- 
tra 

The FUSE spectrum itself is an approximately 
flat continuum, with a mean level of about 3 x 
10~ 14 erg s _1 cm~ 2 A -1 . There are many lines 
of interstellar (IS) H2. No emission lines could 
be seen except those that are definitely associated 
with terrestrial airglow. Most of these lines disap- 
pear in the NE spectrum, with only the strongest 
lines of H 1 and O 1 leaving residual emission peaks 
whose intensities are dramatically reduced com- 
pared to the all photon spectrum. Since the sys- 
temic velocity of PN G 135.9+55.9 is -193 km/sec 
(Tovmassian et al. 2001; Richer et al. 2003), there 
is no possible confusion between airglow emission 
lines and lines from the object. The FUSE NE 
spectrum of PN G 135.9+55.9 is presented in Fig 
1. On the whole, it is impressively featureless, 
particularly when compared with a large subset 
of FUV spectra of planetary nebulae from the 
FUSE archives that we examined. This remark- 
able absence of features in the FUSE spectrum of 
PN G 135.9+55.9 is consistent with the very low 
metal abundances derived from the analysis of the 
PN. 

3.1. Absorption lines from the interstellar 
medium in the FUSE spectrum 

Most of the absorption lines that are observed 
originate in the interstellar matter between us 
and the object. There arc numerous lines of 
molecular and atomic hydrogen. Their study can 
provide interesting constraints upon the distri- 
bution of H and H2 along the line of sight to- 
ward the object (see e.g. McCandliss 2003). How- 
ever, we are mainly concerned with the physics of 
PN G 135.9+55.9 itself in this paper and discuss 
only the impact of these interstellar lines upon our 
analysis of PN G 135.9+55.9. 

Molecular hydrogen is a significant source of 
contamination when studying the spectral features 
of the underlying source. Fortunately, the physics 
of absorption by molecular hydrogen is well un- 
derstood and can be reliably modeled. Absorp- 
tion bands of interstellar molecular hydrogen arise 
from the photo-excitation of molecular hydrogen 
by the background continuum object, causing elec- 



trons in the ground electronic state to be excited 
into higher states of either the Lyman or Werner 
bands. The strength of these absorption lines 
is proportional to the molecular column density 
distribution within the ground electronic state's 
rotational (J") and vibrational {v") 2 energy levels 
along a particular line-of-sight. We used the FUSE 
Simulator Tool to model the FUSE spectrum 
{http://violet.pha.jhu.edu/~gak/fwebsim.html). 
In order to estimate the hydrogen column den- 
sity, we used a flat spectrum of intensity 3 x 
10 -14 erg s _1 cm -2 A -1 and a variety of parame- 
ters for the foreground column densities of atomic 
and molecular hydrogen. Later, we repeated this 
analysis with our derived model spectrum, de- 
scribed in the following sections, and obtained 
exactly the same results. A grid of values of neu- 
tral hydrogen column densities ranging from log 
Ah = 20 to 20.6 cm -2 were investigated. For 
H 2 , the input values were log N(i/', J") = 18 to 
20. The model uses Voigt profiles to calculate 
the optical depth variation with wavelength for 
each line. Line shapes reproducing the linear, 
flat, and square root portions of the curve-of- 
growth are controlled through the specification of 
the Doppler parameter b. The Doppler parameter 
b of line width was adopted equal to 10 km s _1 , 
as provided by default in the simulation tool. 

The modelled interstellar absorption shows a 
significant 0.2A wavelength shift compared to 
the observed spectrum as estimated from cross- 
correlation. This shift is larger than expected for 
interstellar hydrogen absorption, e.g., Bannister et 
al. (2003) did not measure velocities in excess of 
30km/s for interstellar lines from atomic species 
in a sample of white dwarfs with high resolution 
HST UV spectra. Even considering the much 
larger distance to PN G 135.9+55.9 compared to 
that sample, it is unlikely that the entire wave- 
length shift is due to high velocity clouds along 
this line of sight. At least half of the shift in the 
FUSE spectrum of PN G 135.9+55.9 (0.2 A corre- 
sponds to approximately 60 km/sec) is probably 
caused by imperfect centering of the central star in 
the LWRS aperture of FUSE. This is a well-known 
problem that can cause shifts of the wavelength 
scale. We adopted the view that the radial ve- 
locity shift of the H2 features is caused by this 
instrumental effect and corrected the wavelength 
scale for further fitting of the photosphcric line 
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Fig. 1. — The FUSE composite spectrum of PN G 135.9+55.9 for the night only exposure. A few of the 
strongest Hi and Oi A1025.76A airglow lines remain in emission. The other features are absorption lines, 
most of which are due to interstellar matter. Molecular hydrogen lines are marked along the bottom of 
the spectrum: the tick marks for the Lyman series lines are above those for the Werner series. Some other 
prominent interstellar lines of other elements arc marked above the spectrum. The rest positions of Ly/3 and 
Ly7 are marked above the corresponding boxes. 
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profile. 

The H2 column density strongly affects the pro- 
files of the H2 lines, whereas the impact of the 
neutral hydrogen column density is less signifi- 
cant. In Fig. 2 simulated spectra with differ- 
ent hydrogen (Hi, H2) column densities are pre- 
sented. Two spectral regions are shown. The 
first is around A 10 12 A, where the Werner se- 
ries has strong lines, while the other is around 
A 1108 A, with prominent H2 lines of the Lyman 
series. In both regions, even without doing a 
X 2 fit, it is obvious that values of log Nb_ 2 of 
18 and 20 dex are far from fitting the observed 
line profiles and intensities, while logAH 2 = 19 
dex provides a good fit. For neutral hydrogen, 
the choice of the best fit is not so obvious. An- 
other useful web tool allows one to calculate the 
column density depending upon the direction in 
the Galaxy based upon the results of Dickey & 
Lockman (1990, http://heasarc.gsfc.nasa.gov/cgi- 
bin/Tools/w3nh/w3nh.pl). According to this cal- 
culation, the expected neutral hydrogen column 
density is Ah ~ 1.53 x 10 20 cm -2 in the direction 
of PNG 135.9+55.9. This roughly corresponds to 
the modelled value of log Ah = 20.3, which is 
shown in the second panel from the bottom in 
Fig. 2. However a higher value of log 7Vh 2 = 20.6 
seems to be an equally good fit. We suspect that 
the column density toward PN G 135.9+55.9 may 
be slightly higher than 20.3 dex, as follows also 
from the analysis of absorbed flux in the UV. This 
will be discussed in detail below in Sect. 5. 

Besides the numerous H2 lines of the Werner 
and Lyman series, there is a large number of ab- 
sorption lines from O I, N I, N II, N ill and 
Fe 11 that are also due to the intervening interstel- 
lar medium. These lines are commonly encoun- 
tered during FUSE observations (e.g. Froning et 
al. 2001; Wolff et al. 2001) and, most importantly, 
their wavelengths indicate that they are not re- 
lated to our object. The fact that these lines are 
numerous and strong is in agreement with the pre- 
viously deduced large distance to the object and 
its location in the Galactic Halo (Tovmassian et 
al. 2001; Richer et al. 2002, see also Sect. 6). 

3.2. Hydrogen and helium features physi- 
cally related to PN G 135.9+55.9 

The only spectral feature in the FUSE spec- 
trum firmly identified to arise in our object, apart 



from the continuum, is the relatively broad ab- 
sorption from Ly/3 that we attribute to the cen- 
tral star. The only possible emission feature from 
the nebula itself is a small emission peak inside 
the Ly/3 absorption. As mentioned above, we ex- 
pect the lines from the central star to be Doppler 
shifted to the blue by the heliocentric systemic ve- 
locity of the nebula, -193.3 km s _1 (Richer et al. 
2003). The nebular lines may appear to be blue- 
shifted somewhat more, to as much as -205 km 
s , because the nebular emission consists of two 
components with the blue-shifted component be- 
ing the stronger of the two (Richer et al. 2003). 
The portion of the all photon spectrum around 
Ly/3 is presented in Fig. 5 along with the models 
that were fit to it (described below). 

The profile of the absorption feature blueward 
of the vacuum wavelength position of Ly/3, indi- 
cated by a strong airglow line, is extremely com- 
plicated. First, the absorption line is saturated. 
Second, the red wing is contaminated by the air- 
glow emission whose profile itself is complicated 
by a small glitch on the blue wing. Third, there is 
residual emission from the O 1 airglow line further 
to the red. Fourth, there are several strong com- 
ponents of interstellar H2 Lyman lines inside the 
stellar absorption line. Finally, inside the absorp- 
tion feature sits an emission peak that we inter- 
pret as nebular He 11 6-2 A1025.273 emission (see 
below) superposed upon the H 1 Ly/3 absorption 
feature from the central star. The emission line 
that we identify as He 11 has a measured wave- 
length of 1024.57 A. This coincides exactly with 
the position of He II A1025.273 at a radial velocity 
of -205 km s- 1 . 

We corrected the line profile in the all photon 
spectrum, applying the model of interstellar ab- 
sorption with the parameters estimated in the pre- 
vious section (log A H = 20.3 and N R2 = 19). The 
residual spectrum still shows broad photospheric 
Ly/3 absorption as well as emission at the position 
of He 11 line. Although the strength of the absorp- 
tion feature corresponding to Ly/3 depends upon 
the parameters adopted for interstellar hydrogen 
in the model, it is both expected and present in 
the spectrum. We are less confident regarding 
the He 11 emission feature, because it is located 
in a gap between two strong interstellar lines and 
can be easily confused with the residual interven- 
ing continuum. However this supposed emission 
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Fig. 2.— Two portions of all photon FUSE spectrum (LiF 1A and IB channels) of PN G 135.9+55.9 that 
include prominent interstellar lines of H2 are plotted with a dotted line. Ovcrplotted arc modelled spectra 
that assume a flat underlying spectrum of intensity 3 x 10 -14 erg s" 1 cm" 2 A" 1 and different column densities 
of neutral and molecular hydrogen. 



line protrudes more than can be accounted for 
by the models, even if column densities exceed- 
ing logA^H = 20.3 are considered. If this is the 
case, then the upper limit to the flux from the 
He 11 A1025.273 is < 4.4 x 1CT 16 erg cm" 2 s" 1 , ac- 
cording to the model fit that is described in Sect. 
4.2. 

In the optical range, we discovered signatures 
of the central object in the form of absorption 
lines at the wavelengths of hydrogen and helium 
lines. These features were not previously noted 
because our new spectra have higher S/N (by a 
factor of 3) and higher spectral resolution than 
all of our previous spectra. As an experiment, we 
degraded our new CFHT spectra to the resolu- 
tion that we used before, and the absorption fea- 
tures almost disappeared, leaving marginal traces 
only in higher members of the Balmer series. Our 
previous observations also used a wider slit, sub- 
stantially increasing the nebular contribution and 
possibly widening the emission lines. 

As detailed in the Table 1, we obtained seven 
optical spectra of PN G 135.9+55.9 during one 
night. The observations were primarily designed 



to obtain deep spectrophotomctric data in far opti- 
cal UV for study of Ne ions located there and other 
possible faint emission lines from the nebula in the 
integrated spectrum. These will be discussed in 
a subsequent publication concerning the nebular 
characteristics. Here we concentrate on the fea- 
tures due to the central star. Even a simple eye 
inspection of the co-added spectrum immediately 
attracted our attention to the dips around higher 
Balmer series lines, which are clearly not centered 
on the emission lines. An analysis of the individual 
spectra revealed that the center of the absorption 
component varies, indicating large radial velocity 
shifts. 

In principle, our observing technique and equip- 
ment are not ideal for detecting radial velocity 
variations. The use of a grism spectrograph, with 
a nonlinear dispersion, the large range of telescope 
inclination angles spanned by the observations, 
and the variable slit angle combined to scatter the 
absolute zero point of the wavelength scale among 
the spectra by up to 10 A. However, all of our mea- 
surements of absorption components in the central 
star spectrum that are relevant to this study where 
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made relative to corresponding nebular emission 
line components, which effectively eliminates pos- 
sible problems. 

4. Qualitative analysis of the stellar core 
of PNG 135.9+55.9 

4.1. First determinations of the radial ve- 
locity variations 

Initial estimates of the radial velocities (RVs) 
were obtained by de-blending the line profiles with 
the fit routines of the gnuplot program. A Gaus- 
sian profile was assumed for the emission compo- 
nent (from the nebula). A Voigt profile was taken 
for the absorption component, as it is better suited 
to fit a broad white dwarf-like absorption com- 
ponent. We succeeded in separating the compo- 
nents for the lines where we had enough signal, 
enough continuum around the lines, and an ab- 
sorption component sufficiently strong to be re- 
solved. This concerns primarily H7, H<5, and He, 
where de-blending worked fine for all seven spec- 
tra (i.e. for different orbital phases), returning 
consistent parameters of the measured lines from 
one line to another and from spectrum to spec- 
trum. To reduce the number of free parameters, 
we performed a simultaneous fit of all spectra. The 
FWHM of the Gaussian used to fit the emission 
lines and the parameters of the Lorcntzian fit to 
the photospheric absorption were held fixed, be- 
cause these values should remain constant, while 
all other parameters were allowed to vary from 
spectrum to spectrum. 

The imperfections of wavelength calibration 
were eliminated by measuring the radial veloci- 
ties of the absorption components relative to their 
emission counterparts, which we suppose to be 
constant in general, and on such short time scale 
in particular, as is also evident from our previous 
observations (Richer et al. 2002). The correspond- 
ing measurements of the absorption components 
along with epochs of observation are presented in 
the second and sixth columns of Table 1 . Clearly, 
the radial velocity variations of the absorption 
components of these three lines are consistent with 
each other, and vary strongly on a timcscalc of a 
few hours. Before analyzing their time depen- 
dance, i.e. periodic nature, let us consider atmo- 
spheric models of the central star and estimate 
radial velocities from the absorption line fits. 



4.2. Model atmosphere calculations 

The model atmosphere spectra used for our 
analysis were computed with the NLTE code 
PRO 2 developed by Werner (1986). The basic 
assumptions are those of static, plane-parallel at- 
mospheres in hydrostatic and radiative equilib- 
rium. As described in Werner (1986) and Werner 
& Dreizler (1999), the accelerated lambda iter- 
ation (ALI) method is used to solve the set of 
non-linear equations. 

The basic parameters, effective temperature 
and surface gravity, were determined with the grid 
of model atmospheres composed of hydrogen and 
helium calculated by Napiwotzki (1999) for the 
analysis of planetary nebulae central stars. Exten- 
sive model atoms for hydrogen and helium (He I 
and He II) were used (cf. Napiwotzki 1999, for 
details). Line profiles of the hydrogen Balmer 
and Lyman lines and the He II lines were com- 
puted using extended VCS broadening tables (Vi- 
dal, Cooper & Smith 1973) provided by Lcmkc 
(1997) and Schoning & Butler (1989), respectively 

4.3. The gravity of the central star from 
model atmosphere analysis 

It has now become a standard technique to de- 
termine both the temperature and gravity of cen- 
tral stars of planetary nebulae or of hot white 
dwarfs from a simultaneous fit of the Balmer scries 
(e.g. Napiwotzki 1999) or Lyman series (Barstow 
et al. 2001). 

In our FUSE spectrum, only the Ly/3 line has 
a reasonable signal-to-noisc. Higher members of 
the Lyman series fall in the SiC channels, which 
are characterized by lower effective area and hence 
lower S/N. However, the analysis of the Ly/3 line 
is problematic as well, because major parts of the 
line profile are contaminated by geocoronal emis- 
sion or absorption caused by interstellar hydrogen 
molecules. 

The best opportunity for parameter estimation 
is offered by the hydrogen Balmer lines in the op- 
tical spectrum, although these lines are also con- 
taminated by the strong emission coming from 
the nebula. The emission subsides toward higher 
members of the Balmer series, while the decline 
of the photospheric absorption lines is less pro- 
nounced. Thus, it was possible to fit five Balmer 
lines (from H7 to Hg). 
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Table 1: RV measurements from the optical spectra. The heliocentric Julian date of the center of the 
exposures and the RV shifts relative to the nebular emission lines are given b . 



spectrum 


HJD-2,450,000 


RV atm. model fit / 


RV de-blending 


number 




5 lines composite 


H 7 


RS He 


457 


2760.80249 


-204 ± 25 


-95 


-136 -180 


466 


2760.90086 


269 ± 28 


257 


236 365 


467 


2760.92412 


179 ±25 


208 


220 258 


471 


2760.95986 


-53 ± 60 


23 


-7 


472 


2760.98410 


-59 ± 60 


28 


16 -36 


476 


2761.01893 


111 ±28 


178 


138 164 


477 


2761.04169 


246 ± 28 


229 


236 264 



& RV errors estimated from measurements of emission line component does not exceed ±10 km/sec. 



The model atmosphere analysis was done with 
the programme fitsb2, which was originally de- 
veloped to analyze the spectra of double-lined 
white dwarfs from the SPY project Napiwotzki et 
al. (2003, 2004). However, this programme is also 
well-suited for the analysis of the spectra of single- 
lined binaries. fitsb2 performs a simultaneous fit 
of spectra covering different orbital phases, i.e., 
all available information is combined into the pa- 
rameter determination procedure. Fit results arc 
stellar parameters and RVs (a more detailed dis- 
cussion is given in Napiwotzki et al. 2004). Two 
He II lines (4540 A and 4686 A) are visible as well 
and enabled us to estimate the helium abundance. 

We used the nebular emission lines to cor- 
rect the wavelength calibration as described in 
Sect. 4.1. Since the emission components of the 
higher Balmer lines are weak, individual wave- 
length measurements become less accurate. Thus 
we performed an average over all spectra and cal- 
culated wavelength corrections from a simple lin- 
ear fit and applied this for all lines higher than 
H7. Deviations of individual measurements for 
H7, U5, and He, which are most important for the 
RV determination, never exceeded lOkm/s. 

Due to the nebular contamination of the avail- 
able spectra, it was not possible to determine both 
the temperature and gravity from the available 
spectra. Thus, we kept the temperature fixed dur- 
ing the fitting procedure and determined only the 
surface gravity. We considered T c g values in the 
range 80000 K to 150000 K. The results are listed 
in Table 2. Due to potential problems with degra- 
dation of the line profile by orbital smearing (cf. 



discussion below), we restricted the model atmo- 
sphere fit to the three spectra with the largest 
RV shifts (shown in Fig. 3). These were proba- 
bly taken close to quadrature phases, for which the 
rate of change of RV and thus the orbital smearing 
is smallest. 

Below 90000 K and above 120000 K the quality 
of the Balmer line fits degrades significantly. The 
analysis of the nebular spectrum of PN G 135.9±55.9 
indicates a central star temperature of at least 
100,000 K (Tovmassian et al. 2001, Richer et al. 
2002). Taking into account the detection of the 
[Nc V] 3426 A line (Jacoby et al. 2002, Stasihska 
et al in preparation) an even higher temperature is 
required, unless an additional source of EUV pho- 
tons is present in this system. These combined 
constraints are best fulfilled by the solution with 
T eS = 120000 K and log .9 = 5.35. The atmosphere 
model fits to the portions of spectra with largest 
RV shifts are presented in Fig 3. The formal fit 
uncertainty in \ogg for a given value of T c ff is 
0.07 dex. As discussed in Napiwotzki (1999), the 
formal errors computed by the fit procedures are 
usually too optimistic. Thus we made an estimate 
of the external error as described in Napiwotzki 
(1999) and derived an error limit of 0.22 dex. The 
position of PN G 135.9±55.9 in the log g, T e g plane 
is shown in Fig. 4. Although the helium lines arc 
weak and almost completely filled in by the neb- 
ular emission lines, we could derive an estimate 
of the helium abundance (Fig. 3). The fit yielded 
lognHe/«H = —1-0, i.e., close to the value de- 
rived for the nebula. The error of the abundance 
measurement can be estimated from the scatter 
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CFHT 457 CFHT 466 CFHT 477 




i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

20 20 -20 20 -20 20 

AA(A) AA(A) A\(A) 

Fig. 3. — A model atmosphere fit (solid line) to the optical spectra (histogram) taken during the phases with 
largest RV shifts. The dotted lines indicate parts of the spectra which were excluded from the fit. 
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between the individual spectra, which amounts to 
0.2 dex. 




200 100 80 60 40 30 

T e|( / 1000K 



Fig. 4. — Post-AGB evolutionary tracks for differ- 
ent stellar masses from the grid of Blocker (1995) 
and Schonberner (1983) in the (log g, T g) plane. 
The 0.524M Q post-AGB track discussed in the 
text is drawn as dotted line. The connected 
filled symbols indicate the positions of the central 
star of PN G 135.9+55.9 according to the solutions 
ranging from 90000 K to 120000 K. Open symbols 
are central stars of planetary nebulae and other 
post-AGB objects taken from the literature (Napi- 
wotzki 1999) 

The Ly/3 line from the FUSE spectra is sen- 
sitive to gravity as well. However, the total ex- 
posure time of the FUSE spectrum amounts to 
8.5 hours, distributed over 11.7 hours. This is 
much longer than any plausible orbital period 
(Sect 4.4) and it is reasonable to assume that we 
have approximately equal coverage of all orbital 
phases. This means that orbital smearing is the 
dominant broadening mechanism for the observed 
Ly/3 line, even more important than the pressure 
broadening of the photosphcric absorption line. 
We adopted the atmospheric parameters derived 
above (T cff = 120,000,K and log g = 5.35) for the 
fit shown in Fig. 5. 

The photosphcric Ly /? line is contaminated by 
interstellar absorption that we must take it into 
account. In addition to the discussion in Section 
3.1 the column density of the atomic interstellar 
hydrogen can be determined from the red wing of 
the Ly/3 line which is most sensitive to the inter- 
stellar absorption (cf. Fig. 5, note the gap be- 
tween the terrestrial airglow lines of H I and O i). 



Best agreement is achieved with a column density 
Ah = 3 x 10 20 cm -2 . Acceptable fits are derived 
fortherange2xl0 20 cm- 2 < A H i < 5xl0 20 cnr 2 . 
This coincides perfectly with estimates of Ah from 
other wavelength regions (Section 3.1). 

As mentioned above the line widths depend 
upon the assumed amplitude of the orbital mo- 
tion, K\. Wc obtained the best fit for a veloc- 
ity amplitude of 450km/s. Acceptable fits could 
be produced for velocity amplitudes in the range 
from 250km/s to 600km/s. Thus the result from 
the fit of the Ly/3 line is in agreement with the 
orbital solutions discussed in the next section. 

4.4. Investigating the periodicity of the ra- 
dial variations 

As mentioned above, fitsb2 is not only able 
to derive the stellar parameters, but can measure 
RVs as well. We measured RVs for all CFHT 
spectra with temperature and gravity now fixed 
at the value T cS = 120000 K and log g = 5.35 
derived above. The Balmer lines i/7 to H9 and 
the two He II lines were used for the RV deter- 
mination. The results are summarized in Table 1. 
Formal errors provided by the fit procedure are 
relatively small («15km/s). We added to these 
uncertainties an additional lOkm/s to account for 
the remaining uncertainty in the wavelength scale. 
The RV measurements for the conjunction spec- 
tra 471 and 472 are problematic, because the RV 
fits basically rely on the invisibility of the photo- 
spheric Balmer lines. We estimated an RV error 
of 60km/s from the scatter of the results derived 
using different methods. 

The RV variation has a clearly periodic nature. 
However, the sparseness of data does not allow an 
unambiguous determination of the orbital period. 
Fig. 6 shows a "power spectrum" which was pro- 
duced by fitting RV sine-like curves for a range 
of periods. The resulting fit quality, measured 
by % 2 j is plotted in Fig. 6. Since the number of 
data points is low and the measurements are dif- 
ficult due to the strong nebular emission lines, we 
held the system velocity fixed at the value RV = 
-193km/s derived by Richer et al. (2003). The 
best fit quality is obtained for a period of 0^150 
(3^6). The corresponding RV curve is shown in 
the bottom panel of Fig. 7. For this solution, the 
two RV measurements close to conjunction (phase 
0.75) are problematic, because they lie in a "for- 
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Tabic 2: Results of the model atmosphere fits of 
the quadrature spectra 457, 466, and 477 for sev- 
eral values of T e ff. Masses and radii were derived 
from an interpolation in the post-AGB tracks in 
Fig. 4. 



T e fi 

K cm s 



80000 
90000 
100000 
110000 
120000 
130000 
150000 



log g log(n He /"H) M/Mq R/R q 

2 



4.99 
5.05 
5.15 
5.26 
5.35 
5.39 
5.56 



-1.00 



0.62 
0.68 
0.75 
0.83 
0.88 
0.96 
1.02 



0.42 
0.41 

0.38 
0.35 
0.33 
0.33 
0.28 




Fig. 6. — The power spectrum for the FITSB2 
RV measurements from Table 1. The inset gives 
a more detailed view for shorter periods. Periods 
below 0.01 d («15min) can be ruled out, because 
these periods would be much smaller than our ex- 
posure time. 




AA(i) 



Fig. 5. — Model atmosphere fit to the Ly/3 line of 
PNG 135.9+55.9. The histogram is the observed 
all photon spectrum. Only those parts shown with 
a continuous line were used to fit both the stellar 
profile and determine the interstellar absorption. 
The bold, solid line shows the best fit to the stel- 
lar absorption including all of the interstellar con- 
tributions and smearing effects (where K\ is the 
radial velocity of the ionizing star = 450km/s). 
The long dashed lines are the spectra resulting for 
K\ = Okm/s (green) and 900km/s (blue). The ar- 
row indicates the center of the photospheric Ly/3 
line. 



bidden" area of this diagram. These points are 
less problematic for the solutions with shorter pe- 
riods of d 039 (*94) and 0^063 (l h 5), presented in 
the top and middle panels of Fig. 7, respectively. 
Orbital parameters are detailed in Table 3. The 
amplitudes from the best fit solutions for 0^039 d 
and 0^063 arc quite high. We calculated alterna- 
tive ( "minimal" ) solutions by assuming that the 
orbital amplitudes correspond to the maximum 
observed RV shifts (after correction for orbital 
smearing; see below). A solution with a period sig- 
nificantly longer than 4h can be ruled out. Apart 
from the formal x 2 values shown in Fig. 6, simply 
consider the large RV shift between spectra 457 
and 466: 470km/s in just 2^3! 

For the two shortest period solutions found 
(0^063 and 0"?039), note that the exposure times of 
30min extend over a large range of phases. Thus, 
smearing due to the orbital RV variations must be 
taken into account as a major problem. This could 
cast doubt upon whether any absorption line fea- 
tures should be observable at all. We made simple 
simulations of this effect and found that, although 
the line profiles are indeed seriously distorted, line 
profiles similar to those observed can still be pro- 
duced, but that the measured RV values are severe 
underestimates of the real RVs. The difference be- 
tween the real and measured RVs are smallest for 
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P = 0.0391 d 




0.5 0.5 1 1.5 

Phase 



Fig. 7. — Radial velocity curves for three possible 
orbital periods. The error bars in the phase coor- 
dinate correspond to the duration of the individual 
exposures. 

the exposures centered on the quadrature phases 
of largest RV shift. We calculated approximate 
corrections for these phases by measuring the dis- 
placement of the line cores in our simulated spec- 
tra. Results are 124km/s for the 0^039 orbit and 
59km/s for 0^062. The correction for 0"?151 is neg- 
ligible within uncertainties (lOkm/s). We can con- 
clude that the orbital period cannot be longer than 
4 h, but a time series of spectra with much shorter 
exposure times arc needed for an unambiguous pe- 
riod determination. 

It has long been speculated that many plane- 
tary nebulae have binary nuclei. There are models 
that invoke binarity of the central object in order 
to explain the morphology of certain types of plan- 
etary nebulae (Soker 2002). Evolutionary models 
of close binaries require that a fraction of the sys- 
tems going through the common envelope phase 



form a planetary nebula (Yungelson ct al. 1993; 
Iben & Tutukov 1989). However, the number of 
observationally-confirmed binary cores in plane- 
tary nebulae is relatively small. Nonetheless, there 
are reasons to believe that many more of them are 
indeed binaries (De Marco et al. 2003), though 
many probably have relatively long periods. The 
discovery that the central star of PN G 135.9+55.9 
is a close binary is thus very important in this con- 
text. 

5. The extinction of PN G 135.9+55.9 

With our FUSE data, we are able to fur- 
ther address the question of the reddening of 
PNG 135.9 +55.9. In our previous studies, the 
reddening was derived from a comparison of the 
observed hydrogen and helium emission line ra- 
tios with those given by recombination theory. 
Richer et al. (2002) noted that the reddening of 
PNG 135.9+55.9 was modest. The H 7 /H/3 and 
K6/K/3 line ratios gave an E(B - V) of 0.3 - 0.35 
mag (the Ha/H/3 ratio was apparently variable 
and indicated even lower reddening). The He II A 
5412/4686 ratio implied a negative reddening. 

From the FUSE data, we can estimate the ex- 
tinction suffered by the object by fitting the flux 
distribution of PN G 135.9+55.9 from the far UV 
to the near IR. This is shown in Fig. 8. In 
this figure, the observed energy distribution of 
PN G 135.9+55.9 from the far UV to the near IR is 
represented by open squares. The fluxes were mea- 
sured in narrow bands free of any significant lines 
(w 3-4 A in the far UV and w 20 A in the optical). 
The dashed lines are continuum fits to the spec- 
trophotomctric data. The solid line is the model 
of the central star described above for a temper- 
ature of 120 000 K and log g =5.35. The slope 
of the observed UV spectrum (1 000 - 1 200 A) is 
consistent with the model, but shifted downwards 
due to the interstellar extinction. If we correct for 
interstellar reddening using the tables from (Fitz- 
patrick 1999) and assume the standard value i?y 
= 3.1 and E(B — V) = 0.03 mag (Jacoby et al. 
2002), we effectively recover « 0.45 mag of flux, 
but this is not enough to compensate all losses in 
the UV. Increasing E{B — V) to 0.045 mag yields 
a better fit overall, but, at that point, the opti- 
cal fluxes show some excess in comparison to the 
model, while the UV flux is still somewhat short 
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Table 3: Orbital solutions for a 0.55M Q primary and the three possible periods. The set of parameters on 
the left corresponds to the fit results displayed in Fig. 7. For the set of parameters on the right, we assumed 
that the orbital velocity corresponds to the maximum of the measured RV shifts, corrected for the effect of 
orbital smearing. The subscripts 1 and 2 correspond to parameters for the ionizing and companion stars, 
respectively. 



Best fit solution "Minimal" solution 



p 


Vovb 


M% sin i 


a 


Roche 


Roche 


Vovb 


M.% sin i 


a 


^Rochc 


Roche 


(d) 


(km/s) 


(M Q ) 


(Rq) 


(4) 


(4>) 


(km/s) 


(M Q ) 


OR©) 


(Rq) 


OR©) 


0.039 


590 


1.55 


0.63 


0.19 


0.30 


390 


0.73 


0.54 


0.19 


0.22 


0.063 


380 


0.92 


0.77 


0.26 


0.33 


320 


0.69 


0.73 


0.26 


0.29 


0.151 


265 


0.82 


1.35 


0.47 


0.56 













of it. If we also take into account the nebular con- 
tinuum at a temperature of 30 000 K from our best 
fit photoionization models of the nebula (Richer et 
al. 2002, Sect. 8), which contributes primarily in 
the optical range (dash-dot line in Fig. 8), then 
the fit is improved (filled black squares). Much 
better results are obtained using i?v = 2.3 and 
E(B — V) = 0.045 mag (open circles). This could 
mean that the dust towards PN G 135.9+55.9 is 
non-canonical. This may indeed be the case since 
the line of sight toward PN G 135.9+55.9 is not 
necessarily characterized by the canonical value 
Rv = 3.1 (e.g., Barbaro et al. 2001). Moreover, 
the far UV part of Fitzpatrick's reddening curve 
is actually an extrapolation of near UV data into 
the FUSE range. 

Another way to estimate the reddening is to 
use the H I column density derived from our FUSE 
spectrum and convert it into a reddening using the 
relation between E(B — V) and N-r \ of Burstein 
& Heiles (1978). For 7V H i = 3 x 10 20 , we find 
E(B - V) = 0.01 mag. 

These values are somewhat smaller than the 
reddening derived by Richer et al. (2002) from the 
Balmer decrement. Those reddenings were derived 
in the usual way, comparing the observed emission 
line ratios with the theoretical ones and assum- 
ing emission under case B at a nebular tempera- 
ture of 30 000 K. However, no account was taken 
of the underlying stellar absorption. When this 
effect is properly taken into account (Stasihska et 
al. in preparation), the reddening derived from 
the Balmer decrement is compatible with the val- 
ues obtained from our FUSE spectrum. 

Thus we confirm that the reddening, and hence 



extinction, of PN G 135.9+55.9 is indeed low. In 
the following, we adopt a value of E(B — V) = 0.04 
mag, but none of our conclusions would be signif- 
icantly affected should the real value be slightly 
different. 

6. The distance and age of PN G 135.9+55.9 

Having an estimate of the stellar gravity and 
a measurement of the apparent magnitude of the 
central star, we can estimate the distance d to the 
object using the expression: 

d= l.ll(M*F v 10°- 4mv /ff) - 5 

where d is in kiloparsecs, M+ is the stellar mass 
in solar units, i*V is the model stellar flux in the 
V band in 10 8 erg cm~ 2 s _1 A -1 , and my is 
the observed stellar magnitude corrected for neb- 
ular contribution and for extinction. This method 
to evaluate the distances to PNe has been dis- 
cussed extensively by Napiwotzki (2001) and there 
are strong arguments that such spectroscopic dis- 
tances are reliable (within the uncertainties). To 
evaluate the distance to our object, we take the 
observed magnitude my = 17.9 mag from Richer 
et al.(2002), assume a reddening of E(B — V) = 
0.04 mag, and apply a correction of 0.15 mag for 
the nebular continuum contribution, as derived 
from photoionization models described in Richer 
et al. (2002). The stellar flux for the model with 
T off = 120000 K and log 5 = 5.35 is 7.07 • 10 8 
erg cm~ 2 s _1 A -1 . Adopting a stellar mass of 
0.55M Q , we find a distance to PN G 135.9+55.9 
of 17.8 kpc. Had we adopted a stellar mass of 
0.88M Q , as indicated by Fig. 4 and listed in Table 
2, we would have found a distance of 22.6 kpc. 
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Fig. 8. — The flux distribution of PN G 135.9+55.9 from far UV to near IR. The open squares are the 
observed fluxes. The dashed lines are continuum fits to the observed data. Fluxes corrected for interstellar 
extinction are represented by filled squares (Ry =3.1, E{B — V) = 0.045 mag). The dot-dash line is the 
expected contribution from the nebula (see text) . The open circles are the extinction-corrected data (Ry = 
2.3) after subtraction of the nebular component. The central star model with 120 000 K and \ogg = 5.35 is 
presented as a solid line. The model is normalized tqghe raw flux observed at 4 000 A. 



Assuming that the stellar temperature is only 
90000 K and taking the values of \ogg and M* 
from Table 2, we find d = 24.2 kpc, while assum- 
ing a stellar temperature of 150000 K, we find d 
= 20.6 kpc. The largest contribution to the error 
budget is the uncertainty in log g for a fixed T e fj, 
which amounts to an uncertainty of 5 kpc in the 
distance. 

In any case, our present estimates of the dis- 
tance are consistent with PN G 135.9+55.9 being 
located in the Galactic halo. Note that the values 
of the distance obtained in this way are roughly 
consistent with the ones obtained previously (Tov- 
massian et al. 2001; Richer et al. 2002), which were 
based upon nebular properties (flux and dimen- 
sions), whereas the present determination relies 
only on the fit of the stellar spectrum. This indi- 
cates that we are progressively getting a coherent 
picture of the nebula. 

On the other hand, the discovery of the com- 
pact binary nature of the nucleus of PN G 135.9+55.9 
introduces new challenges in interpretation of the 
evolutionary path of this unique object. For fur- 
ther discussion, a rough estimate of the age of the 
planetary nebula is useful. The expansion time, 
^cxp, of the planetary nebula is generally computed 
as iexp = -Rout/^oxp, where i? ut is the outer radius 
of the nebula and i; oxp is its expansion velocity. 
This is of course a simplistic approach, since it 
is known that the expansion velocity should vary 
during the course of the nebular evolution (see 
e.g. the dynamical models of Villaver, Manchado 
& Garcfa-Scgura 2002) and that the epoch of the 
ejection of the nebula does not necessarily coin- 
cide with the evolution of the central star off the 
AGB. In the case of PN G 135.9+55.9 the expan- 
sion time measures the time since a final mass loss 
episode triggered by common envelope interaction 
of the binary stellar core (Sect. 7). The transition 
phase from the AGB to a hot post- AGB was cer- 
tainly strongly affected by the binary evolution. 
Nevertheless, an estimate of the nebula expansion 
time provides a valuable constraint, which is es- 
pecially important in the context of close binary 
stellar evolution. We adopt the nebular size de- 
termined by Richer et al. (2002) and an expansion 
velocity of 30 km s _1 as found by Richer et al. 
(2003). For a distance of 17.8 kpc, we find a neb- 
ular expansion time of 15800 yr. Using this rough 
age estimate, the position of the object in the 
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Fig. 9. — Post- AGB evolutionary tracks for dif- 
ferent stellar masses interpolated from grid of 
Blocker (1995); Schonberner (1983) in the (T cff , 
tcxp) plane. Tracks are plotted with solid/dotted 
lines for the parts before/after maximum temper- 
ature. The position of PN G 135.9+55.9 is marked 
for the estimated values of T e ff = 120 000 K and 
t cxp = 16 000yr. 

(T e ff, t exp ) diagram (Figure 9) can be compared 
with that of model stellar tracks of various masses 
(Blocker 1995; Schonberner 1983). The position 
of PNG 135.9+55.9 in Fig. 4 indicates that the 
comparison must be made with the part of the 
tracks that lie prior to the temperature maximum 
(solid lines in Fig. 9). We find that the position of 
PN G 135.9+55.9 is consistent with a mass of the 
ionizing star around 0.57M Q . 

7. Discussion: The status of the stellar 
core of PN G 135.9+55.9 

7.1. The ionizing star 

As described above, we derived the surface 
gravity of PN G 135.9+55.9 from fits of the Balmer 
lines of selected optical spectra. Because of nebu- 
lar contamination of the optical spectra, we cannot 
determine both the temperature and gravity, so we 
adopted T e g = 120000 K based upon the analysis 
of the nebula (Stasihska et al. in preparation). 
The resulting gravity is logg = 5.35. From a 
comparison with evolutionary tracks of H-burning 
central stars of planetary nebulae (Fig. 4) we esti- 
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mated a mass of 0.88M Q . This value is at variance 
with the expected properties of a Population II 
central star. Three lines of evidence point to- 
wards a Population II nature of PN G 135.9+55.9: 
the large height above the Galactic plane (z = 
14.8kpc for d = 17.8 kpc), the high radial velocity 
(RV = — 193km/s), and the very low metal abun- 
dances derived from the analysis of the planetary 
nebula and the far-UV spectrum of the central star 
(Stasihska et al., in preparation). 

Known post-AGB stars of Population II have 
low masses (« O.55M ; e.g., Moehler et al. 1998), 
which is explained by the small mass of their pro- 
genitors. One could argue that problems caused 
by the orbital motion and the contamination by 
the nebula emission lines cause a systematic shift 
of the fit results. However, from a comparison of 
observed and synthetic spectra, we conclude that a 
gravity much higher than that given by our fits can 
be ruled out, because the resulting lines would be 
much broader than observed. Our estimate of the 
effective temperature is based upon the analysis of 
the nebular spectrum. If an additional source of 
hard ionizing photons were present in this system, 
(e.g. resulting from interaction between the com- 
ponents of the binary system), the nebular spec- 
trum would actually overestimate the temperature 
of the ionizing star. Lower temperatures would 
result in lower mass estimates (Table 2) in better 
agreement with the expected mass of a Popula- 
tion II central star. 

Due to the obvious Population II nature of 
PN G 135.9+55.9, in the subsequent discussion we 
shall adopt a mass of 0.55M Q for the ionizing star. 
A low mass is also supported by the relatively large 
kinematic age of the nebula. If the mass of the 
ionizing star were really as high as the estimates 
derived from Fig. 4, e.g., 0.9 M©, the time interval 
between the tip of the AGB and the hottest point 
of the post-AGB evolution would be of order of 
only 100 years, much shorter than the estimated 
nebular expansion age of 16000 years. The ob- 
served nebular age is in much better agreement 
with a mass of « 0.57Mq (or lower), for which 
transition times arc «20000yrs (or longer). Note 
that the common envelope phase very likely ab- 
breviated the post-AGB evolution. 

How might we explain the apparent discrepancy 
between the observed and expected parameters of 
the ionizing star? The progenitor almost certainly 



went through a common envelope phase, which 
drove it out of thermal equilibrium. In this case, 
the relations used to derive the stellar mass from 
log g and T c g do not apply. However, relaxation 
happens on the thermal time scale of the post- 
AGB envelope, which is given by 



where G is the gravitational constant, M the mass 
of the central star, M env the mass of the stellar 
envelope, R the stellar radius, and L the stel- 
lar luminosity. For the derived stellar parameters 
and an adopted mass for the post-AGB envelope 
of Mcnv = 10 _2 M©, the resulting thermal time 
scale amounts to only «30 years. Since the real 
envelope mass is probably smaller (cf. Fig. 3 of 
Schonberncr 1983), this can be considered an up- 
per limit. The short relaxation time makes it un- 
likely that we observe the star precisely in this 
stage. 

A more plausible explanation is that the post- 
AGB star experienced a late thermal pulse (shell 
helium flash) after emerging from the common en- 
velope phase. Stars suffering a late thermal pulse 
evolve back to cooler temperatures and then repli- 
cate the previous post-AGB evolution (Iben et al. 
1983; Schonberner 1983; Blocker 1995). The evo- 
lutionary speed and luminosity during this sec- 
ond phase can be quite different from the ones 
in the first phase. As an illustration, we plotted 
in Fig. 4 the track computed by Blocker (1995) 
for a 0.52M© post-AGB star, which experiences 
two late thermal pulses. The resulting track spans 
a wide range of gravities and luminosities, whose 
details depend upon the precise mass and thermal 
pulse phase. Most importantly, the time scales 
for this kind of evolution are thousands to tens of 
thousands of years, making it much more likely to 
catch a star in one of these phases. 

7.2. The companion star 

The short orbital period and large scmiampli- 
tude of the RV that we detect imply that the other 
component of the binary system should be a rather 
massive star. The mass function can be used to 
derive lower limits to the mass of the unseen com- 
panion, if the mass of the visible star is known 
(Table 3). In the case of the 3^6 solution we de- 
rive a minimum companion mass of 0.82M© for a 



17 



O.55M ionizing star (1.O1M for a O.88M ioniz- 
ing star). The best fit orbits for the shorter period 
solutions indicate even higher companion masses. 
However, the "minimal" solutions yield lower lim- 
its of «0.70Mq. Thus the companion must be a 
relatively massive star. 

Note that the system is physically small. The 
separation between both stars amounts to only 
1.35i?0 for the 3^6 solution and the 0.55Mq ion- 
izing star (Table 3) . The Roche lobe of the com- 
panion has a size of 0.56i?Q. A metal-poor 0.6Af Q 
main sequence star has a radius of 0.55-R© (Baraffe 
et al. 1997), almost filling its Roche lobe if the 
period is 3.6 h and overfilling it for the other solu- 
tions. However, the orbital solutions require more 
massive companions that have even larger radii 
and would therefore not fit into their Roche lobes, 
so a main sequence companion appears unlikely. 
Further evidence against a main sequence compan- 
ion comes from the the binary system BE UMa. It 
consists of a hot post-AGB star with parameters 
similar to that of PN G 135.9+55.9 (Liebert et al. 
1995) and a K dwarf companion. However, this 
system has a much wider orbit (P = 2*?29) than 
PN G 135.9+55.9, yet it still shows a spectacular 
reflection effect producing a forest of RV variable 
emission lines from the heated surface of the com- 
panion (Ferguson & James 1994). The absence of 
these lines in our spectra is another strong argu- 
ment in favor of a compact companion, i.e., a white 
dwarf or a neutron star. A photometric monitor- 
ing campaign to search for variability should be 
able to provide final evidence against (or for) a 
main sequence companion. Note that the ioniz- 
ing star fills a significant fraction of its Roche lobe 
as well (Tables 2 and 3). The T cff = 120000 K, 
logg = 5.35 solution corresponds to a radius of 
O.26i?0 compared to a Roche lobe of 0.47i?© in the 
3.6 h orbit. According to the temperature-radius 
relation from the 0.546Af Q track of Schonberner 
(1983) the ionizing star completely filled its Roche 
lobe when it had a temperature of 50000 K. This 
represents the minimum temperature of the ion- 
izing star for the stellar core of PN G 135.9+55.9 
to become a detached system. Model predictions 
indicate that the central star fills or overfills its 
Roche lobe for the 1.6 h and 0.94 h solutions, re- 
spectively. 

In the standard picture of the common enve- 
lope phase the companion spirals in until enough 



orbital energy is released to eject the envelope 
(Wcbbink 1984). For a system with a compan- 
ion as massive as in the case of PNG 135.9+55.9, 
one would expect that only a moderate spiral in 
would be necessary, producing a relatively wide 
system. We used standard assumptions for the 
common envelope phase (formula A. 14 of Nele- 
mans et al. 2001, with oceA = 2) to estimate 
the initial separation of the PNG 135.9+55.9 sys- 
tem before the final common envelope. For the 
3^6 orbit we calculated an initial separation of 
only 3.6i?0 (initial/final mass: O.85M Q /O.55M0, 
companion mass: 0.8M@). This is an impossible 
value, because this would indicate that interaction 
started shortly after the main sequence phase, pre- 
venting the star from becoming the central star of 
a planetary nebula. Thus it is difficult to explain 
the formation of the PN G 135.9+55.9 system with 
standard assumptions. On the other hand this 
means that PN G 135.9+55.9 could provide an im- 
portant test of our understanding of common en- 
velope evolution. 

If the companion is a compact star, the binary 
will merge in less than 1 Gyr due to the loss of 
energy and angular momentum caused by gravi- 
tational wave radiation. Our estimates from the 
mass function show that the companion should 
be rather massive and it is quite possible that 
the sum of the masses of this system exceeds the 
Chandrasckhar limit for white dwarfs (1.4M©). 
Such systems are proposed as progenitors of SN la 
(Webbink 1984; Iben & Tutukov 1989). The SPY 
project (Supernovae la Progenitor surveY) Napi- 
wotzki et al. (2001, 2003) performed a radial ve- 
locity survey of white dwarfs for potential SN la 
progenitors. To date, only one possible candidate 
system was detected (Napiwotzki et al. 2003), il- 
lustrating how rare these systems are. If it can be 
shown that the core of PN G 135.9+55.9 is mas- 
sive enough to qualify as SNIa progenitor, this 
would provide the bonus that additional informa- 
tion about the formation of such systems is avail- 
able from the chemistry and morphology of the 
planetary nebula. 

8. Conclusions 

We have presented new optical and far-UV 
spectra of the central star of the very metal-poor 
planetary nebula PN G 135.9+55.9. The photo- 



18 



spheric Balmcr lines visible in the optical spec- 
tra enabled us to determine the surface gravity 
from a model atmosphere analysis. For a tem- 
perature of T e ff = 120000 K (as indicated by the 
nebular emission lines), we derived a gravity of 
logg = 5.35 and a photosphcric helium abundance 
lognne/^H = —1-0. Our model atmosphere spec- 
trum also fits the UV Ly/3 line profile observed in 
the FUSE spectrum satisfactorily given that the 
Ly/3 line is contaminated by interstellar lines and 
suffers from severe smearing effects due to the or- 
bital RV variations. 

Comparing the above values of T e g and logg 
with those from evolutionary tracks in the logg- 
T c g plane yields a mass of 0.88M Q for the ionizing 
star, if one adopts standard H-burning post-AGB 
stellar tracks. This result is in contrast with the 
Population II nature of this object, indicated by 
its metallicity (Tovmassian et al. 2001, Richer 
et al 2002, Jacoby et al. 2002, Stasihska et al., 
in preparation), location in the galactic halo, and 
its large radial velocity (Richer et al. 2003). We 
argue that the real mass of the ionizing star of 
PN G 135.9+55.9 is close to the typical masses of 
Population II post-AGB stars: w 0.55M Q . Since 
the stellar core of PNG 135.9+55.9 is a close bi- 
nary the formation of the planetary nebula is likely 
to have happened during a phase of common en- 
velope evolution. 

A comparison of a model spectrum calculated 
for appropriate stellar parameters and the ob- 
served FUSE and optical fluxes indicates a red- 
dening of Eb-v — 0.04 mag. This relatively 
small value is in agreement with the expecta- 
tions for an object at high galactic latitudes 
(b = +56°). This low reddening is also in agree- 
ment with the estimates the interstellar atomic 
and molecular hydrogen column densities towards 
PNG 135.9+55.9 based upon the FUSE observa- 
tions. Allowing for extinction, we obtain a dis- 
tance for PN G 135.9+55.9 from a comparison of 
its measured flux in the V band and model atmo- 
sphere fluxes. The result is «18kpc, which places 
PN G 135.9+55.9 15kpc above the galactic plane. 
The kincmatical age of the nebula is «16000 years, 
in rough agreement with expectations if the cen- 
tral star mass is of order 0.55M Q . This expansion 
age argues against a more massive central star, 
e.g., 0.9Mq, since the expected post-AGB age of 
such a star is only 100 years. 



The optical spectra reveal large radial veloc- 
ity variations of the stellar absorption lines, obvi- 
ously caused by orbital movement in a close binary 
system. Since the observations were not designed 
for the measurement of a binary orbit, the phase 
coverage and the wavelength calibration are not 
optimal. However, we show that the orbital pe- 
riod must be less than 4 hours (our three best 
solutions are 0.94, 1.5, and 3^6 hours) and the 
minimum (half-)amplitudc of the radial velocity 
amounts to 250km/s. The companion remains 
invisible in the spectra during all orbital phases. 
Our preliminary orbital solutions indicate a com- 
panion mass of « 0.9Mq. A solid lower limit is 
O.6M0. It is very unlikely that the companion is a 
main sequence star, because we would then expect 
a strong reflection effect producing very strong 
emission lines from the irradiated hemisphere of 
the main sequence star. That these are not ob- 
served in the optical spectra indicates that the 
companion is most probably a white dwarf or a 
neutron star. 

Due to the short orbital period (four hours or 
less) the system will merge in less than one Gyr, 
if the companion is a compact star. If a compan- 
ion mass exceeding 0.85M Q can be confirmed, the 
total mass of the stellar binary core would exceed 
the Chandrasekhar limit and this system would 
qualify as a potential progenitor of a type la su- 
pernova. New observations, with shorter exposure 
times and longer time coverage, will be needed in 
order to establish a definitive orbital period and 
other parameters for the binary system. 
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